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Abstract 

\ In this paper we suggest an iterative algorithm to compute automatically the 

■ scattering matrix elements of any given effective lagrangian, T. By exploiting the 

relation between T and the connected Green function generator, Z, we provide a 
formula which does not require the use of the Feynman graphs and it is suitable 
to implement a numerical routine. By means of this algorithm we have built a 
relatively simple and fast fortran code which we have used to calculate, at the tree 
level, the rate of four fermion production at LEP II (finding a very good agreement 
with previous calculation) with and without the emission of one observable photon. 



Introduction 

The computation of the matrix element of physical processes is often a difficult task. 
With the increasing collider energies the complexity of the required computations has 
grown considerably. When the number of final particles is high, even to evaluate the 
corresponding tree level Feynman diagrams becomes hard and the final formula is often 
an intricated function of several variables, inadequate for practical use. 

For this reason, many authors made a lot of work to simplify these calculations. Very 
good achievements have been obtained: the helicity amplitude technique Q shows how to 
manipulate the Feynman rules to compute by hand some extremely difficult multi-particle 
processes which otherwise would be impossible to handle, color decomposition, recursion 
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relations || and string inspired |§ rearrangement of the diagrammatic expansion have 
allowed a drastic simplification in the number of graphs to be considered in some QCD 
processes, programs of symbolic manipulation provide a very useful tool to manage and 
evaluate the Feynman graphs Q and finally a completely automatic treatment of the 
Feynman rules is available ||. 

Here we pursue a different approach. By exploiting the relation among the one particle 
irreducible Green Functions generator F with the connected Green Functions generator 
Z we are able to provide an explicit formula for Z (for a given T) as a series expansion 
and we show that after the truncation of the series to a proper number of steps the exact 
scattering amplitude is recovered. We use this formula to implement a FORTRAN code 
and we present the computation of the rates for the processes e + e~ — *> /i + z/ At r _ z/ T and 
e + e~ — > fi + v fl T~i? T 'y for some center of mass energies accessible at LEP E. 



An iterative algorithm for the scattering matrix 

If the one particle irreducible Green Functions generator, T, of a theory is known the 
computation of the S-matrix requires the evaluation of the legendre transform, Z, of T M 

Z(ji, ...,j n ) = -T((p u ...,<p n ) +j k (x)(p k (x) (1) 

where are the classical fields defined as the solution of 

J k = f- (2) 

and the play the role of classical sources. 

For the sake of clarity, in the following we will stick to a specific example, namely we 
discuss the process 

<K<?i) + <f>(<h) -> <K<? 3 ) + 0(g 4 ) (3) 

in a A0 3 theory. As a matter of convenience we use the convention that all the particles 
are incoming and therefore the momenta of the outcoming particle take an overall minus 
sign. From now on we define 

Pi =Qi P2 = <?2 P3 = -93 Pa = -<74 (4) 
The lagrangianP] £ of the model, is 

C = J d'x{- l -<pd^ + ^) 

" P rd 4 q \p 2 5(p + q)<P(p)<P(q) + ^ / ,f P , A f t \ d A k 8(p + q + k)(f){p)(f){q)(f){k). 



(2tt) 4 * T ^ ,;rwrw/ 6 J (27r) 4 (2vr) z 

(5) 



1 In the following example we take the most simple lagrangian, but the generalization to an arbitrary 
effective lagrangian (see after) which also includes fermions and higher order operators (e.g. radiative 
corrections terms calculated by hand) is almost straightforward. 
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Hence the Z functional for the theory is, at the tree level, the Legendre transform of (g) 



Z = -£ + J d 4 x (f>(x)J{x) (6) 

In order to calculate the scattering matrix element for the process (§) we specify the 
source J(x) to 



J(x)=J>e^ (7) 
1=1 

where a>i is a coefficient which might carry lorentz and internal symmetry indices. The 
required matrix element is obtained after differentiation of Z with respect to the a; in 
the limit a\ — > and after having properly cut the external legs. With the source J(x) of 
eq. (|7|) the functional Z becomes 

Z(ai, a 2 , a 3 , a 4 ) = —C + / d A x ^ a m e iPm ' x 4>(x) (8) 

J m=l,4 

where is such that it minimizes (or more generally it picks up a stationary point) 
the right hand side of eq. (||) for any fixed set of a m , i.e. it is of the form 

<j>(x) = b m e- iP ™ x (9) 

with 

p m = c l mPl 4 = o, +i. (io) 

In fact, even if the the fulfillment of the stationary condition for an arbitrary choice of the 
a,,- would require c l m to take any integer value, we are interested in the Green functions 
which are recovered only after taking the limit a,j — > 0. In such a limit, these are the only 
momenta which survive and which contribute to the final result^. In our specific case we 
have fourteen Pj (and bj) 

Pb=P\+P2] Pe=Pi+P3; 

Pl3 = Pi +P3 +Pi\ Pl4 = P2 + P3 + Pi\ 

Replacing <p{x) in eq. (||) we obtain, in momentum space, 

Z(a 1 , a 2 , a 3 , a 4 ) = ~Y1 ~ \ bjbib m D jlm + ^ (H) 

j,l j,l,m m,l *m 



2 Even if theoretically one should take only the real part of such J{x), in practice, for our pragmatic 
aim this is not significant. 

3 After, when the explicit iterative solution for the stationary condition will be given, it will be manifest 
that other momenta will show up only with higher order aj . 



3 



where the and A; m take into account the momentum conservation, i.e. 





= A 


if Pi + Pj + P k = 




= 


if P l + P j + P k ^0 


A m l 


= P 2 

m 


if P m + P ; = 


A m i 


= 


if P m + Pi^0. 



and the bj variables in eq. ( |TT| ) satisfy the minimum equations 

Auk + ^D ijk bjbk £ = 1,4 

Atf&j + ^Dijkbjh £ > 4 (13) 



and 



— 

dbi 



Here an important point should be noticed: we began to study the problem of finding the 
functional Z of eq. @ which is, formally, a problem with an infinite number of degrees 
of freedom. Using eq. (0) with the constraint given in eq. (10) we have reduced it to the 



solution of eqs. fljjD where the number of degrees of freedom is finite. 

We now solve iteratively eq. (p~3|) by replacing the bj with the sum bj = Y^Lo bj,r where 
the bj tr are iteratively defined as functions of the bj jS , with s < r, as follows: 

1 = b ii0 £ = 1,4 

= b ii0 i>A (14) 



where we set dj = p 2 to obtain the truncated Green function. At the second step the 
same eq. fll3|) gives 

1 

" L 

2 



bj,i — —-AjlDikmbkfibmfl (15) 



and finally as last step 



bj2 — —-AADikmibkflbm^ + &fc,l& m ,o) (16) 



To converge to the exact solution of eq. (|13D one should continue the series of the & Jir up 
to r = oo, however, for the purpose of computing the scattering amplitude, the terms 
we evaluated in (|l4],|l^,|n|) are sufficient. In fact the bj^ series can be regarded as an 
expansion in the coupling A and we have indeed truncated the series to the order A 2 . 

We now substitute the obtained bj >r in eq. (pTT|) . The result coincides with the required 
scattering amplitude 

2 2 2 2 <9 4 zT(ai, a 2 , 03, Q4) 
A PlP2P3Pi - Pl p 2 p 3 p 4 tim daida2da3daA 
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i 14 i 14 

_ o E E b j,rbi, s Aji — p XI E Djimbj,rbi >s b m ,t 

Z s+r=2j,l=l D s+r+t=l j,l,m=l 

4 14 

i=i j=i 

x2, 1 1 1 



(jPl + P2) 2 (Pi + P3) 2 (Pi + P±f 



(17) 
(18) 



One must notice that eq. fll7|) contains some redundacy and if one is wishing to use it to 
build a numerical routine he would immediately face the problem to handle quantities of 
the type 0/0. It is however possible to use eqs. (p~3|) to simplify eq. ([17]) (and simultane- 
ously overcoming the problem of dealing with 0/0 expressions). Multiplying both sides of 
eq. (|13D by bi and summing over the index i, one obtains, from eq. (17), 



2 2 2 2 d 4 ^(ai, a 2 , 03, 04) _ 1 v y- D-, b- b, b * (19) 

1 z * s,r,t j,(,m=l 

where the indices j,l,m are defined according to the conventional ordering given above. 
Notice that the evaluation of eq. ( |T§| ) does not require the knowledge of the 6 J)2 of eq. ( |TE| ) 
which are non zero only if j > 10. Due to this fact one does not need to evaluate the 
inverse of the propagator for on shell four momenta which occurs in eq. (16). 



For an arbitrary lagrangian the previous iterations can be generalized in the following 
way Q 

2 6 

b h = Xj j = l,n 

b" = j > n n = number of external particle 
h a -Irr^n^ 75 V K 1 b 5 

U j,m — Q ■?>* MM U k,r U l,s 

r+s=m— 1 



Z s+r=n-2 s+r+t=n-3 



(20) 



where greek symbols denote lorentz or internal symmetry indices, <ft denotes a generic field 
of the theory (a sum over distinct particles is understood), x? stand for a properly chosen 
source term (for instance , \ = 1 for a scalar, \ a = u a (p) for a Dirac spinor with four 
components a, etc.), II is the appropriate propagator, II is the inverse of II and O is 
the generalization of the Djm of the eq. ([3~3|) which can depend on the spin, the momenta 
and any other internal quantum number. In eq. fl20D we have included only trilinear 

4 In eqs. (f20| ) we have implicitly assumed that all the fields <jf are real. For charged fields, one can, 
obviously, use the more conventional (and practical) complex notation. We did not do it in order to avoid 
a tedious digression about the combinatorial factors. 
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interactions; notice, however, that any operator with a higher dimension can always be 
reduced to a set of trilinear and bilinear ones using a proper set of auxiliary fields. 

A final remark is in order here. In eqs. ( p0[) one has to take properly into account 
particle statistic. Formally this can be achieved introducing a set of Grassman variables 
6j, €j€k + = 0, and setting u(p) — > eu(p) for the fermion sources where now u{p) is a 
vector of ordinary numbers. In practice this means that each term in the sums of eqs. (P0| ) 
enters with a relative sign depending on the order of the bj. 

Finally one should keep in mind that a simplification closely parallel to the one ob- 
tained in eq. ( |I9|) do always occur. 

The advantage of the iterative eqs. (p0|) is that they can now be easily implemented 
in a fortran code. 



ALPHA, an ALgorithm to compute PHysical Amplitudes 



We have implemented a fortran code, ALPHA, based on the solution of eqs. (PUD- We have 
first studied a purely scalar theory which allows to check also processes with an high num- 
ber of particles. We have then added fermions with a Yukawa type interaction and verified 
that the fermionic statistic is properly accounted for. Finally the complete electroweak 
Lagrangian, in the unitary gauge, has been added (at present the strong interactions are 
absent). A very important test for the code is provided by the gauge properties of the 
theory. The good behavior of the scattering amplitudes, calculated by ALPHA, with the 
increasing center of mass energies shows that delicate gauge cancellations do occur. 

We have then compared our results with known one and reproduced the cross sections 
for many processes with two particles in the final state and the rates for the processes 
e + e~ — > WWZ, e + e~ — > ZZZ with a good accuracy. 

We now want to add some technical details about the code. To stress again the simplic- 
ity of eqs. (p0|) we remark that the bulk of the code consists of less than 1000 FORTRAN 
executable statements. On top of this one has to add several peripheral subroutines, for 
example subroutines to read the input, to return the Dirac 7 matrix, to perform four vec- 
tor scalar products, etc. All these peripheral subroutines have two common features: they 
are completely trivial and they can be checked one by one in a completely independent 
way and therefore they do not increase the complexity of debugging. 

Another important feature, which is again linked to the simplicity of eqs. (p0[), is the 



efficiency achieved by the code. As a practical demonstration of this fact we show in table 
[I] the CPU performance for various processes. It is clearly seen that the increment in 
CPU time is significantly smaller than the increment in the number of Feynman graphs 
and in their complexity. Furthermore ALPHA is only a first prototype of the algorithm 
and we believe that a large room for optimization is left. 



Emission of a hard photon in four fermion production at LEP H 

In this section we will use our technique to compute the cross section of four fermion 
production from e + e~ collision at the LEP I center of mass energy and when an additional 
hard photon is found in the final state. 
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process 


cpu time (seconds) 


e + e~ — > n+v^r^Ur 


128 


e + e~ -> / u + z/ M r"z/ r 7 


170 


e + e~ — > yU + z/ M r _ z/ T 77 


395 


e + e~ — > e + z/ e e"z/ e 


180 


e + e~ — > e + z/ e e"i/ e 7 


323 


e + e~ — > e + z/ e e"z/ e 77 


877 



Table 1: Required CPU time for the computation (single precision) of 10000 matrix 
elements for various processes. All the Higgs couplings are set to zero. The times are 
given in seconds and the calculations have been performed with a DIGITAL machine 
ALPHA 3000/600 with 64M of memory. 



One important goal of the LEP/SLC collider program has been achieved by the precise 
measurements. After the first phase which has revealed an important and enlightening 
agreement for most |7], H] of the experimental data and the Standard Model predictions 
|J, the measurement of H^-boson mass at LEP I will complete the pattern of precision 
tests. The importance of this goal has been pointed out by several authors The 
most promising way to measure the charged boson mass at LEP II is to make a precise 
measurement of the cross section near above the threshold of H^-pair production. The 
two body phase space is proportional to the velocity (3 of the two final H^-boson; thus the 
total cross section near the threshold essentially behaves as (3. The strong dependence 
of this parameter from the boson mass means that the measurement of the total cross 
section can be seen as a precise measurement of the H^-mass. 

Experimentally the iy-bosons can be detected only through their decay products. As 
a first approximation the rate for each process is given by the total VT-pair cross section 
times the branching ratios of W into a pair of fermions, but the existence of a finite width 
and of other non-resonant diagrams introduce sizable corrections. Due to the experimental 
accuracy which is aimed to at LEP II the full computation of the process 

e + + — > 4 fermions (21) 

is required^ |TT|. In the first column of table (Q) we show the results of the program 
ALPHA for the process 

e + + e~ — ► /i + + z/ M + t~ + v T (22) 

near above the PU-pairs production threshold, and in the second column the set of the 
relevant input constants is shown. We have applied a cut over the invariant mass of 
Ufj, and u T of 3 GEV to eliminate the contribution from r pair production with subsequent 
r decay. 



5 However, when proper kinematical cuts are applied, there are semianalytical approximations which 
reproduce the correct result very accurately 
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y/S in GeV 


cross sect, in pb 


input parameter 


value 


150 
157 
160 
163 
170 
180 
200 


0.00594 (3) 
0.00937 (4) 
0.04201 (14) 
0.09113 (26) 
0.17019 (52) 
0.21166 (47) 
0.22891 (55) 


M z 
sin 2 (6 w ) 
<Xqe D (2M w ) 
M w 
T w 
T z 


91.1888 GeV 
0.23103 
128.07 
80.23 GeV 
2.03367 GeV 
2.4974 GeV 
1.7 GeV 



Table 2: W-pair production rates (in pb) and the input constants. 

The value of the sin 2 {6 w ) and of the a~^ ED (2M w ) given in the table enter into the 
cross section through the tree level formula of the coupling describing the fermion and 
gauge boson SM interactions. Since the W-boson lifetime is finite, the cross section is 
only smoothly suppressed below the S = ^M^ threshold; it increases up to the 200 GeV 
region and after decreases as 1/5. The results are in agreement within the montecarlo 



error, with the existing computations for the process (22|). 

Beyond the born level contribution, the QED virtual and real corrections are the most 
significant. Soft or collinear photons cannot be seen by real experimental detector, thus 
only a cross section that includes the emission of such photons can be compared with the 



experimental results [11] 



In the following we will not deal with the problem of infrared and collinear QED 
correction to the process (p2"|) and we will present the results of the study of the exclusive 
process 

e + + e" — ► /x + + + t~ + v T + 7- (23) 

The number of Feynman diagrams increases considerably and it is very difficult to 
handle the calculation with the ordinary methods. The matrix element is a function of 
eleven variables and also the numerical integration of such a big number of variables is 
not easy. 

To achieve the numerical accuracy for each total cross section shown in table (f|), a 
set of about 1.7 millions of events has been collected. To perform the integral over the 
kinematical variables we have used the package VEGA ]l2j to properly take into account 
the peaking behavior of scattering matrix element. 

We have introduced two cut-off to keep the infrared divergences under control and to 
provide experimentally testable results: 

- E 7 > lGev is the minimum photon energy required; 

- 6 > 200 mrad (cos (6) < 0.98) is the minimum angle allowed between the photon 
and each (both initial and final) charged fermion in the process. 

The table (|3]) shows the total cross section at various initial energies. If an integrated 
luminosity of about 500 pb~ x is collected, few events will pass the above cuts with a 
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in GeV 


cross sect, in fb 


number of events ( 500pb r ) 


150 


0.464 (2) 


0.23 


157 


1.300 (4) 


0.65 


160 


2.874 (12) 


1.44 


163 


5.724 (17) 


2.35 


170 


14.795 (56) 


7.4 


180 


20.638 (52) 


10.3 


200 


24.070 (62) 


12.0 



Table 3: r + v T /i 7 production rates: E y > 1 GeV and 9 > 200 mrad (see the text). 

differential distribution as in figures (1) and (2). It is clear from figure (1) that there is a 
forbidden region (i.e. too low rate) and a region where some events will show up. 

Conclusion 

In this paper we have suggested an iterative algorithm which we have found very simple 
and effective to compute scattering amplitudes. 

Differently from other approaches it never requires the explicit use of Feynman graphs 
since the truncated connected Green functions are directly derived by performing a nu- 
merical Legendre transform of the effective Lagrangian. We also have shown that the 
algorithm can be successfully used to implement a fortran code, that can calculate very 
cumbersome processes in a relatively short time. To provide an example of the perfor- 
mances of our code we have computed the cross sections for some relevant processes of 
physical interest during the LEP I phase, namely the processes e + e~ — > h + v^t~u t and 
e + e~ — >• /z + z/ M r~z/ r 7. 
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Figure Captions 

Figure 1: differential cross section (fb/BIN) of e + e~ — > /i + z/ M r~z/ T 7 as function of 0.42 < 
cos{6) < 0.98 and 1 GeV < £ 7 < 19.5 Gev (for E cm = 200 GeV). 

Figure 2: differential cross section (fb/BIN) of e + e~ — > fi + v tl T~v T / -f as function of cos(#) 
(for £ cm = 200 GeV). £ 7 > 1 GeV (see the text). 
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